MCNP6 Fission Cross Section Calculations 
at Intermediate and High Energies 
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MCNP6 has been Validated and Verified (V&V) against intermediate- and high-energy fission- 
cross-section experimental data. An error in the calculation of fission cross sections of 181 Ta and a 
few nearby target nuclei by the CEM03.03 event generator in MCNP6 and al "bug" in the calcu- 
lation of fission cross sections with the GENXS option of MCNP6 while using the LAQGSM03.03 
event generator were detected during our V&V work. After fixing both problems, we find that 
MCNP6 using CEM03.03 and LAQGSM03.03 calculates fission cross sections in good agreement 
with available experimental data for reactions induced by nucleons, pions, and photons on both 
subactinide and actinide nuclei at incident energies from several tens of MeV to about 1 TeV. 
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MCNP6 [U is used in various applications involving fis- 
sion reactions at low, but also at intermediate and high 
energies. It is critical that it describes such reactions 
as well as possible, therefore it is often validated and 
verified against available experimental data and calcu- 
lations by other models (see, e.g., [l|, and refrences 
therein). However, for nuclear reactions involving fission, 
generally the transport codes arc only tested for how well 
they describe the fission-fragment yields, and the emis- 
sion of particles from them (spectra and multiplicities, 
both "prompt" and "delayed"). Calculation of fission 
cross sections, af, themselves by the transport codes is 
usually not tested, as it is assumed that data libraries 
used at energies below 150 MeV must provide reliable 
<jf, as should the eventgenerators used at higher ener- 
gies. Our recent work [3J shows that this assumption is 
too optimistic, as some event generators have problems 
describing well some fission cross sections. In such cases, 
all the characteristics of the corresponding fission reac- 
tions calculated by the transport codes, like the yields of 
fission fragments, spectra and multiplicities of neutrons 
and other particles, both "prompt" and "delayed" would 
also not be reliable. To assess this situation, we test how 
MCNP6 calculates <jf using the Cascade- Exciton Model 
(CEM) and the Los Alamos version of the Quark-Gluon 
String Model (LAQGSM) as realized in the event gen- 
erators CEM03.03 and LAQGSM03.03, respectively, for 
reactions induced by nucleons, pions, and photons at in- 
cident energies from several tens of MeV to ~ 1 TeV. 
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A comprehensive description of the CEM03.03 and 
LAQGSM03.03 event generators can be found in Ref. 
Q Therefore, we present only a brief summary of how 
CEM03.03 and LAQGSM03.03 calculate fission cross sec- 
tions, needed to better understand our results. 



II. a f CALCULATION BY CEM AND LAQGSM 

CEM03.03 and LAQGSM03.03 calculate the fission 
cross sections (and the yields of fission products and evap- 
oration of particles from them) using a modification of the 
Generalized-Evaporation-Model code GEM2 by Furihata 
[5j. GEM2, in its turn, is based on a modification of the 
older RAL fission-model code by Atchison (see and 
references therein). For fissioning nuclei with Zj < 88, 
GEM2 uses the original Atchison calculation of the neu- 
tron emission width T 
the fission probability as 



and fission width F f to estimate 



Pf = 



f 



1 



r/ + r„ i + r„/r y 



(i) 



Atchison uses [f| the Weisskopf and Ewing statistical 
model 0- All details on the RAL and GEM2 codes and 
all formulas used by them to calculate cr/ can be found 
in Refs. [f| @. Let us mention here only that in the 
case of subactinide nuclei, the main parameter that de- 
termines the fission cross sections calculated by GEM2 
is the level-density parameter in the fission channel, a/ 
(or more exactly, the ratio a//a n ), where a n is the level- 
density parameter for neutron evaporation. GEM2 uses 
for it the follwing approximation: 



a f 



1.08926 + 0.01098(x - 31.08551) 5 



(2) 
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with x = Z 2 /A. 
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FIG. 1. Prokofiev systematics (open circles) and exper- 
imental proton-induced fission cross sections of 181 Ta (sym- 
bols, see detailed references in [j|) compared with our old 
MCNP6 calculations (black dashed lines) using the CEM03.03 
event generator before we fixed the error, with the corrected 
CEM03.03 results (green line), and with calculations by the 
updated MCNP6 using the corrected CEM03.03 event gener- 
ator (red dashed line), as indicated. 



Furihata obtained Eq. (2) as the best approximation 
for the ratio a//a„ while calculating 07 with GEM2 used 
after the Bertini IntraNuclear Cascade (INC) model @, 
without taking into account possible preequilibrium emis- 
sion of particles after the INC. The INC's of CEM and 
LAQGSM are different from the Bertini version Q, and 
both models account for preequilibrium emission of par- 
ticles. Hence, for any particular reaction, the mean mass, 

< A >, and charge, < Z >, numbers of compound nuclei 
that may fission, as well as their mean excitation energy, 

< E >, calculated by CEM and LAQGSM differ from the 
ones provided by the Bertini INC. As a result, the best 
approximation provided by Eq. (2) for the Bertini INC 
is not the best for CEM and LAQGSM. In Ref. j|, we 
performed our own fitting of the ratio a t j a n which works 
the best for CEM and LAQGSM, using as "experimen- 
tal" proton-induced 07 the systematics by Prokofiev (loj . 
After this adjustment for proton-induced reactions, CEM 
and LAQGSM were tested on reactions induced by neu- 
trons, photons, and pions, and were found to calculate 
well 07 for these types of reactions (see details in @). 



III. RESULTS 

Rowing from [§] that both CEM and LAQGSM de- 
scribe well a/ for various reactions, we did not check 
until very recently how well MCNP6 and its precursor, 
MCNPX, using CEM and LAQGSM, calculated fission 
cross sections. The first attempt to calculate 0/ with 
MCNP6 using CEM03.03 for p + 181 Ta found some very 
bad results at proton energies above ~ 200 MeV (see 
the black dashed line in Fig. 1). We verified that these 



MCNP6 results were not a problem of the incorporation 
of CEM03.03 into MCNP6, but were a result of an error 
inadvertently introduced into CEM03.03 in 2005 during 
a major upgrade of the code [ll| . 
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FIG. 2. Prokofiev systematics [l(| (open circles) and exper- 
imental proton-induced fission cross sections of 181 Ta (sym- 
bols, see detailed references in 0]) compared with our results 
by the updated MCNP6 using the LAQGSM03.03 event gen- 
erator with the GENXS option (red dashed lines) and with 
calculations by LAQGSM03.03 used as a stand alone code 
(green lines), as indicated. For comparison, wrong initial re- 
sults by an older version of MCNP6 (called "Beta 1"), before 
the "bug" in the calculation of the fission cross section with 
the GENXS option while using LAQGSM03.03 was fixed, are 
shown as well with several magenta crosses. 

To fix this problem, we have refitted in Q the values 
of at/a n in GEM2 for 181 Ta. Our new CEM03.03 results 
for 181 Ta are shown in Fig. 1 with a green line. We have 
replaced in MCNP6 the initial defective CEM03.03 mod- 
ule with the corrected version with the correct values of 
a f /a n in GEM2. Results from the updated MCNP6 arc 
also shown in Fig. 1 by a red dashed line. 

Having discovered this 2005 error and knowing how 
it affects the CEM results, we can understand now why 
in a recent work 12J it was found that CEM03.02 (which 
has practically the same physics as the version CEM03.03 
used here) provided such a poor agreement with the mea- 
sured yields of the nuclides produced in proton interac- 
tions with 181 Ta and nearby target nuclei for energies 
above 250 MeV. 

Unfortunately, we also discovered Q that 07 printed in 
the output files of MCNP6 while using the GENXS op- 
tion with LAQGSM03.03 were not correct; we present an 
example of such erroneous results in Fig. 2, shown in the 
plot with magenta crosses. MCNP6 provides correct re- 
sults for fission-fragment yields and for particle spectra in 
its output file. However, an unobserved error in counting 
the number of fissioning nuclei in the GENXS portion of 
MCNP6 when using the LAQGSM03.03 event generator 
was present in the initial Beta 1, version of MCNP6 
We have fixed that error. The current version of MCNP6 
is free of that error and provides in its output files correct 
values for <jf (see the red dashed line in Fig. 2). 
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FIG. 3. Examples of 07 for reactions induced by p, n, 7, and n~ . Detailed references to measured data (symbols) can by found 
in [j§|. Dfferences in the absolute values of 07 calculated with CEM03.03 used as a stand-alone code and with MCNP6 using 
CEM are related to different total reaction cross sections calculated by CEM and used by MCNP6 for these reactions. 



After fixing both these problems, we find Q that 
MCNP6 using CEM03.03 and LAQGSM03.03 calculates 
<7/ in good agreement with available data for reactions in- 
duced by p, n, 7, and tt~ on both subactinide and actinidc 
nuclei, at incident energies from several tens of MeV up 
to <~ 1 TeV (see an example in Fig. 3 and more details 
inRef. 0). 

IV. CONCLUSIONS 

MCNP6 has been validated and verified against 
intermediate- and high-energy fission-cross-section exper- 



imental data. An error in the calculation of fission cross 
sections of 181 Ta and other nearby target nuclei by the 
CEM03.03 event generator of MCNP6 and a "bug" in the 
calculation of fission cross sections with the GENXS op- 
tion of MCNP6 while using the LAQGSM03.03 event gen- 
erator were detected during our current V&V work. After 
fixing both these problems, we find that MCNP6 using 
CEM03.03 and LAQGSM03.03 event generators calcu- 
lates a/ in a good agreement with available experimental 
data for reactions induced by nucleons, pions, and pho- 
tons on both subactinide and actinide nuclei (from 165 Ho 
to 239 Pu) at incident energies from several tens of MeV 
to about 1 TeV. 
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